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In the present paper we construct a self-consistent theory, interpreting the observational properties of RBS1774. It is well known 
that the distribution function of relativistic particles is one-dimensional at the pulsar surface. However, cyclotron instability causes 
an appearance of transverse momenta of relativistic electrons, which as a result, start to radiate in the synchrotron regime. We study 
the process of the quasi-linear diffusion developed by means of the cyclotron instability on the basis of the Vlasov's kinetic equation. 
This mechanism provides generation of measured optical and X-ray emission on the light cylinder lengthscales. A different approach 
of the synchrotron theory is considered, giving the spectral energy distribution that is in a good agreement with the XMM-Newton 
observational data. We also provide the possible explanation of the spectral feature at 0.7keV, in the framework of the model. 
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1. Introduction 

RBS1774 (1RXS J214303.7+065419) has been the most re- 
cent XDIN (Xray di m isolated neutron star) to be found 
(IZampieri et al. L l200ll) . Its X-ray spectrum is well reproduced 
by an absorbed blackbody with a temperature kT ~ lOOeV and 
with a total column density of Hh ~ 3 • 10 20 cra~ 2 . Application of 
more sophisticated, and physically motivated models for the sur- 
face emission (atmospher i c mod els) result in w orse agreement 
with th e data dZane et al. L I2005I) . According to ISchwope et al. I 
(2009), a fit to the X-ray spectra extracted from RGS spectro- 
graphs onboard XMM-Newton yields that the best result is ob- 
tained when the two-temperature blackbody model is used. But 
the same model applied to the X-ray spectra extracted from three 
EPIC detectors does not improve the fit compared to the simple 
blackbody model. However, the formation of a non-uniform dis- 
tribution of the surface temperature is more likely artificial and 
needs to be examined by convincing theory. 

Alternatively, the observational properties of RBS1774 can 
be explained i n the framework of the pla sma emi ssion model first 
develo ped by [Machabel i & Usovl (1 19791) and lLominadze et al. I 
, d!983l) . According to these works, in the electron-positron 
plasma of a pulsar magnetosphere the waves excited by the cy- 
clotron resonance interact with particles, leading to the appear- 
ance of pitch angles, which obviously causes synchrotron ra- 
diation. We suppose that the X-ray emission from this object 
is generated by the synchrotron mechanism. According to the 
stand a rd theory of the s ynchrotron emission dBekefi & Barrett L 
Il977t iGinzburgl. 1 198 lh the typical synchrotron spectrum is a 
power-law, when the present model suggests different spectral 
distribution. The main reason for this is that we take into account 
the mechanism of creation of the pitch angles, consequently re- 
stricting their values. Contrary to this, in the standard theory of 
the synchrotron radiation, it is assumed that along the line of 
sight the magnetic field is chaotic, leading to the broad inter- 
val (from to n) of the pitch angles. The present model gives 
successful fit for the observed X-ray spectrum, when the orig- 
inally excited cyclotron modes enter the same domain as the 



measured optical emission of RBS1774. We suppose that the 
observed spectral feature at 0.7keV in the X-ray spectrum of 
RBS1774 is caused by wave damping process developed near 
the light cylinder due to the cyclotron instability. 

In this paper, we describe the emission model (Sec. 2), de- 
rive theoretical X-ray spectrum of RBS1774 and fit with XMM- 
Newton observations (Sec. 3), explain the possible nature of the 
observed spectral feature at ~ 0.7keV (Sec. 4), and discuss our 
results (Sec. 5). 

2. Emission model 

The distribution function of relativistic particles is one dimen- 
sional at the pulsar surface, because any transverse momenta 
(Pi.) °f relativistic electrons are lost in a very short time(< 
10" 20 s) via synchrotron emission in very strong magnetic fields. 
For typical pulsars the plasma consists of the following compo- 
nents: the bulk of plasma with an average Lorentz-factor y p 
10 2 , a tail on the distribution function with y, 10 s , and the pri- 
mary beam with jb — 10 7 (see Fig . [TJ . However, plasma with an 
anisotropic distribution function becomes unstable, which can 
lead to a wave excitation in the pulsar magnetosphere. The gen- 
eration of waves is possible during the further motion of the rela- 
tivistic particles along the dipolar magne tic field lines if the con - 
dition of cyclotron resonance is fulfilled (Kazbegi et aT71 ll991h : 



a>-k,V„ 



, Cl>B 

k x u x H 

Jr 



0. 



(1) 



where u x = cV^jrjpojB is the drift velocity of the particles due to 
curvature of the field lines, p is the radius of curvature of the field 
lines and cl>b — eB/mc is the cyclotron frequency. During the 
generation of waves by resonant particles, one als o has a simulta- 
neous feedback of these waves on the electrons dVedenov et all 
fl96lt) . This mechanism is described by quasi-linear diffusion, 
leading to the diffusion of particles as along as across the mag- 
netic field lines. Therefore, resonant electrons acquire transverse 
momenta (pitch angles) and, as a result, start to radiate through 
the synchrotron mechanism. 
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Fig. 1. Distribution function of a one-dimensional plasma in the 
pulsar magnetosphere. Left corresponds to secondary particles, 
right to the primary beam. 

The kinetic equation for the distri bution function of the res- 
onant particles can be written as dMachabeli & UsovT 1 19791; 
iMachabelieuO 120021) : 

where G is the force responsible for conserving the adiabatic 
invariant p\/B(r) =const, F- is the radiation deceleration force 
produced by synchrotron emission, and g|| is the reaction force 
of the curvature radiation. They can be written in the form: 



mc 

G± = Jrlft, 

P 



G\\ = —Jrty , 



F x = -a s ^(l+y 2 r i^), F,| 



2 , 2 

-a s y r ip , 



(3) 
(4) 

Gii = -cxcyt, (5) 

where a s = 2e 2 a> 2 B /3c 2 , a c = 2e 2 /3p 2 and if/ m p±/p\\ <s 1 is the 
pitch angle. 

Now let us compare the transverse components of the forces 
G and F. If we consider the case yip » 1 we will have : 



G± 3m i c 6 y I r 



F ± 2e 4 B 2 p\R s ) rV 2 



1 



(6) 



where B s is the magnetic field at the star surface, R s is the star 
radius and r is the distance from the pulsar. For the typical pa- 
rameter values of pulsars |G_J <K \F ± \. Then taking into account 
that d/dif/ » d/dy the equation for th e diffusion across the mag- 
netic field can be written in the form (Chkheidz e et al. L 12010) 



at p ± dp ± v ' p x op ± \ dp ± j 

where 

D ± , ± = -6\E k \ 2 (8) 

is the diffusion coefficient and \Ek? is the density of electric en- 
ergy in the waves. 

The transversal quasi-linear diffusion increases the pitch- 
angle, whereas force F resists this process, leading to the sta- 
tionary state {df/dt = 0). Then the solution of Eq. (7) is jy ± . 



f x = Cexp 
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(9) 



To evaluate p ±0 , we use the quantity 

2 co 

where a> is the frequency of original waves, excited during the 
cyclotron resonance and can be estimated fromEq. (1) as follows 
co ~ cosldy,-. Consequently, we will get 
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The mean value of the pitch-angle fa as p± /p\\ - 10 -3 (i. e. the 
assumption y r i//o » 1 done at the beginning of our computations 
proves to be true). As a result of the appearance of the pitch 
angles, the synchrotron emission is generated. 

3. X-ray spectrum 

Let us consider the synchrotron emission of the set of elec- 
trons. The number of emitting particles in the elementary dV 
volume is p ±f dp ± dp\\dVdQ. T , with momenta from the intervals 
[p±,P± + dp±] and [pn,/>n + dp\\], and with the velocities that 
lie inside the solid angle d£l T near the direction of r. If we 
write the parallel distribution function of the emitting particles as 
J p±fd p± =./J | (pn), then t he emission flux of the set of electrons 
will be (Ginzburgl [T98ll) 



Ief\\(P\\)dp\\dVdQ. T , 



(12) 



where l e is the Stokes parameter, which is additive in this case, 
as the observed synchrotron radiation wavelength A is much less 
than the value of rT 1 ^ - the average distance between particles, 
where n is the density of plasma c omponent electro ns. The inte- 
gral (fT2l is easily reduced to (see lGinzbur^l d 19811) ) 



f/^iiwf r 

J £ m YJel 



dp\\ 



(13) 



K 5/ i(z)dz 

e/e m J 

Here e m ~ 5 ■ lO^fii/o^keV is the photon energy of the maxi- 
mum of synchrotron spectrum of a single electron and Ks/^z) is 
a Macdonald function. After substituting the mean value of the 
pitch-angle in the above expression for e,„, we get 



£,n -5-10- 



1/9 5 7 9 \ 1/4 

r 1 ' 2 ( 3m c y h i 



32e 6 P 3 j 



(14) 



For the primary beam electrons with the Lorentz factor yi, ~ 10 7 
the emitted photon energy e m ~ O.lkeV comes in the energy 
domain of the observed X-ray emission of RBS 1774. Thus we 
suppose that the measured X-ray spectrum is the result of the 
synchrotron emission of primary beam electrons (the resonance 
occurs on the right slope of the distribution function of beam 
electrons (see Fig. |TJ), switched on as the result of acquirement 
of pitch angles by particles during the quasi-linear stage of the 
cyclotron instability. 

To find the synchrotron flux in our case, we need to know the 
one-dimensional distribution function of the emitting particles/j|. 
Let us multiply both sides of Eq. (2) on p ± and integrate it over 
p ± . Using Eqs. (3), (4 ), (5) and the following expressions for the 
diffusion coefficients (Chkh eidze et al. 11201 Oh 



(15) 
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And also taking into account that the distribution function van- 
ishes at the boundaries of integration, Eq. (2) reduces to 



d A = JL 

dt dp\\ 



.2^1/2 ^-Lo 



2 4 

Pa* + a c y 



4mc 2 y 



(16) 



Let us estimate the contribution of different terms on the right- 
hand side of Eq. (16). The estimations show that the first term 
is much bigger than two other terms. Consequently, for the 
primary -beam electrons instead of Eq. (16), one gets 



d A 8 l ° s 2 A 

dt d Pn \m 2 c 2 7i l i 2P±Jl r 



dp 

Considering the quasi-stationary case we find 
1 



p\ !2 m 



(17) 



(18) 



For yifj 10 10 , a magnetic field inhomogeneity does not 
affect the process of wave excitation. The equation that de- 
scribes the cyclotron no ise level, in this case, has the form 
dLominadze et al. Lfl983b 



d\E k \ 2 

dt 
where 



= 2T c \E k \% 



2 2 

n e 

E c — — : — f\\(Pres)i 



(19) 



(20) 



is the growth rate of the instability. Here k\\ can be found from 
the resonance condition ( 1 ) 



Km ss . 

c5y res 

Combining Eqs. (17) and (19) one finds 



d_ 

dt 



d P\\ 



„V2 



a = 



.6 ,,3 \ 



4 e z a u B r p 

3 7t 5 C 5 U 2 j 



1/4 



which reduces to 



m 

I/ 2 



= const. 



(21) 



(22) 



(23) 



(24) 



Taking into account that for the initial moment the major 
contribution of the lefthand side of the Eq. (24) comes from/j| o , 
the corresponding expression writes as 



dPW 



„ i i 2 



-f\\o- 



(25) 



The distribution function/is proportional to n ~ 1/r 3 , then one 
should neglect /j| in comparison with /f| n . Consequently, the above 
equation reduces to 



d P\\ 



,V2 



+4 = o. 



(26) 



As we can see the function Ekipw) drastically depends on the 
form of the initial distribution of the primary beam electrons. 
Here we assume that the initial energy distribution in the beam 
has a Gaussian shape 



fbo ~ 



exp 



(y - jb) 



(27) 



where yj ^ 10 - is the half width of the distribution function 
and rib = B/Pce is the density of p rimary beam electrons, equa l 
to the Goldreich- Julian density (iGoldreich & Julian! [l969). 
Since jj <K ji,, this distribution is very close to <5-function. 
Consequently, the electron distribution can be taken as monoen- 
ergetic. 

In this case for the energy density of the waves we get 
\E k \ 2 K Pll . (28) 
The effective value of the pitch angle depends on \E k \ 2 as follows 



1 

2lo e 



a 2 3 , ? V /4 

; T\Ek\ 

P\ J P 



(29) 



According to our emission model, the observed radiation 
comes from a region near the light cylinder radius, where the 
magn etic field lines are pract ically straight and parallel to each 
other (lOsmanov et al. L [2009), therefore, electrons with ijj « ipo 
efficiently emit in the observer's direction. 

Using expression (18), (28) and (29), and replacing the inte- 
gration variable p\\ by x — e/e m , from Eq. (13) we will get 



F f oc e 



K 5/3 (z)dz 



dx. 



(30) 



The energy of the beam electrons vary in a small interval. In 
this case the integral (30) can be approximately expressed by the 
following function 



F e oc e 03 exp(-e/e m ). 



(31) 



We performed a spectral analysis by fitting the model spec- 
trum (Eq. (31)) absorbed by cold interstellar matter, with the 
combined data extracted of the three EPIC X-ray cameras of 
the XMM-Newton Telescope. The resulting x 2 — 1 -63 and the 
amount of interstellar matter «// = (3.36±0.2)xl0 20 cm~ 2 , which 
appears to be close to the tota l Galactic absorption in the source 
direction {n„ = 5 • \0 20 cmT 2 [Pickev & Lockman I (1 19901) ). The 
spectral feature at ~ 0.7keV that is mostly described as an ab- 
sorption edge or line (IZane et al. 1 120051: ISchwope et all [2009) 
is also evident in our case from inspection of Fig. [2] We find that 
adding an absorption edge at 0.7keV improves the fit, leading 
to a reduced x 2 — 1 -50. The best-fitting energy of the edge is 
Eedge = 0.679keV, and the optical depth is T e d ge = 0.20 (see 
Fig. [3). The fitting results are listed in Table 1. We suppose 
that existence of the absorption feature in spectra of RBS 1774 
is caused by wave damping at photon energies ~ 0.7keV, which 
takes place near the light cylinder. 

4. Possible nature of the spectral feature 

During the farther motion in the pulsar magnetosphere, the X- 
ray emission of RBS 1774 that is generated on the light cylin- 
der lengfhscales, might come in the cyclotron damping range 
dKhechinashvil i & Meli kidze 111 9971) : 



U)B 

oj — k.. V, - k x u x = 0. 

Jr 



(32) 
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The condition for the development of the cyclotron instability 
may be easily derived for the small angles of propagation with 
respect to the magnetic field. Representing the dispersion of the 
waves as 



co = kc, 



(33) 



and neglecting the drift term, the resonance condition (32) may 
then be written as 



1 

lf r + 2 



LOB 

ary r 



(34) 



where 6 a if/ is the angle between the wave vector and the mag- 
netic field. Taking into account that ip^ » l/2y^ one finds from 
Eq. (34) the frequency of damped waves 



2u) h 



m = 



2 ' 



(35) 



If we assume that the resonant particles are the primary beam 
electrons, then the estimation shows that on the light cylinder 
lengthscales eo = (/z/27r)<Jo - 0.7keV. 

5. Discussion 

According to the generally accepted point of view, the X-ray 
spectrum of RBS 1774 is purely thermal and is best repre- 
sented by a Planckian shape. A fit with a pure blackbody com- 
ponent absorbed by c old interstellar matter gives x 1 - 1-81 
dSchwope et al. L l2009t) . Including a Gaussian absorption line at 
~ 0.7keV (as the largest discrepansies between model and data 
are around 0.7keV) improves the fit x 2 — 1 -50 (parameters are 
listed in Tab. 1). However, the nature of this spectral feature is 
not fully clarified as yet. The most likely interpretation is that it 
is due to proton cyclotron resonan ce, which implies ultras trong 
magn etic field of B cyc ~ 10 14 G dZane et al. L 120051; iRea et al. L 
2007). Although, the required strong magnetic field is inconstis- 

tent with timing measurements giving Bn p = 3.2 • 10 19 VPP 
2 • 10 13 G ( Ka plan & van KerkwiikLl2009h . 



We are not about to reject the existing thermal emission mod- 
els, but in present paper we propose an alternative explanation of 
the observed X-ray spectrum of RBS 1774. It is supposed that the 
emission of this source is generated by the synchrotron mecha- 
nism. The distribution function of relativistic particles is one di- 
mensional at the pulsar surface, but plasma with an anisotropic 
distribution function is unstable which can lead to wave excita- 
tion. The main mechanism of wave generation in plasmas of the 
pulsar magnetosphere is the cyclotron instability, which devel- 
ops on the light cylinder lengthscales. During the quasi-linear 
stage of the instability, a diffusion of particles arises along and 
across the magnetic field lines. Therefore, plasma particles ac- 
quire transverse momenta and, as a result, the synchrotron mech- 
anism is switched on. If the resonant particles are the primary 
beam electrons with yi, 10 7 their synchrotron emission en- 
ter the same energy domain as the measured X-ray spectrum of 
RBS 1774. 

We construct a self-consistent theory interpreting the obser- 
vations of RBS 1774. Differently from th e standard theory of the 
synchrotron emission (Ginzburg , 1981), which only provides a 
power-law spectrum with the spectral index less than 1, our ap- 
proach gives the possibility to obtain different spectral energy 
distributions. In the standard theory of the synchrotron emis- 
sion, it is supposed that the observed radiation is collected from 
a large spacial region in various parts of which, the magnetic 
field is oriented randomly. Thus, it is supposed that along the 
line of sight the magnetic field directions are chaotic and when 
finding emission flux, Eq. (13) is averaged over all directions of 
the magnetic field (which means integration over i// varying from 
to 7r). In our case the emission comes from a region of the pul- 
sar magnetosphere where the magnetic field lines are practically 
straight and parallel to each other. And differently from standard 
theory, we take into account the mechanism of creation of the 
pitch angles. Thus we obtain a certain distribution function of 
the emitting particles from their perpendicular momenta (see Eq. 
(10)), which restricts the possible values of the pitch angles. In 
the framework of the model we obtain the following theoretical 
X-ray spectrum of RBS 1774 F £ oc e 03 exp(-e/e m ) and perform 
a spectral analysis by fitting data from the three EPIC detectors 
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Table 1. The model parameters of RBS 1774 for combined fits to EPIC-pn and EPIC-M OS in the energy i nterva l 0.2 -1.5 keV(The 
fitting results with a pure blackbody model absorbed by cold interstellar matter are from Schwop e et al. I (120091) 1 



TV/ToHpI 

1V1UUC1 


(10 2( W 2 ) 




(eV) 




K1 bb 
(eV) 


^edge/line 

(eV) 


&~ line 

(eV) 


^"edge/line 




plasma 


3 36 +0 - 20 


7.0 


±0.2 












1.63(311) 


plasma*edge 


3.30+°;j 2 


6.9 


±0.1 






679 + _\l 




20 +0 03 

"•■^"-0.03 


1.50(309) 


bbody 


1.85*" 






103.5 


±0.8 








1.81(311) 


bbody*gabs 


1.84**° 






105.1 


±0.9 


131% 


97+I6 
z/ -4 


6 5 +12 


1.50(308) 



simultaneously. The fit with a model spectrum absorbed by cold 
interstellar matter yields^ 2 = 1.63 (see parameters in Tab. 1). 

During the farther motion in the pulsar magnetosphere, the 
X-ray emission of RBS 1774 comes in the cyclotron damping 
range (see Eq. (32)). If we assume that damping happens on the 
left slope of the distribution function of primary beam electrons 
(see Fig. [TJ, then the photon energy of damped waves will be 
eo = (h/2n)(L>o - (hl2n)2cL>Bljb^ 2 — 0.7keV. Taking into ac- 
count the shape of the distribution function of beam electrons, 
we interpret the large residuals around ~ 0.7keV (see Fig. [2j as 
an absorption edge. Including an absorption edge improves the 
fit leading to a reduced^ 2 = 1.50. The best-fitting energy of the 
edge is E ei i ge - 0.679keV, and the optical depth is r e d ge = 0.20 
(see Tab.l). However, adding an absorption edge to the model 
spectrum does not produce a sta tistically signifi c ant im prove- 
ment of the fitting. According to Schwope et al. I d2009) if one 
uses the RGS X-ray spectra of RBS 1774 in place of EPIC spec- 
tra, the resulting^- 2 is changed just marginally when a Gaussian 
absorption line is included at ~ 0.7keV. Thus, we conclude that 
the nature of the feature at 0.7keV is uncertain and might be 
related to calibration uncertainties of the CCDs and the RGS at 
those very soft X-ray energies. The same can be told about a fea- 
ture at ~ 0.3keV (the large residuals around 0.3keV are evident 
from inspection of Fig. [2] and [3}. A feature of possible simi- 
lar nature was detected in EPIC-pn spectra of the much brighter 
prototypical object RX J 1856.4-3754 a nd classified as remaining 
calibration problem bv lHaberTI (120071) . Consequently, more data 
are necessary to finally prove or disprove the existence of those 
features. 

The frequency of the original waves, excited during the cy- 
clotron resonance can be estimated from Eq. (1) as follows 
v » 2jTO)B/6jb ~ 10 I4 Hz. As we can see the frequency of cy- 
clotron modes comes in the same domain as the measured op- 
tical e mission of RBS 1774 (IZane et al. L l2008t ISchwope et all 
2009). 
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